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The models of star formation function and of dissipation of turbulent energy of interstellar medium are 
^ ' proposed. In star formation model the feedback of supernovae is taken into account. It is shown that 
Q , hierarchical scenario of galaxy formation with proposed models of star formation and dissipation is able to 
Q ' explain the observable star formation pause in the Galaxy. 
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B ' Introduction 



In models of evolution of galaxies where the direct gas-dynamical computation is used, it is difficult to 

^ , reach the resolution in space less than 10 — 100 pc and resolution in time less than, say, 10^ yr due to 

■ computational cost. We can use only the phenomenological approach to describe the interstellar medium 

I— i" (ISM) on scale 100 pc and below. This means that we have to use an averaged description of interstellar 

^ . medium. 

Q-i' In the present paper we consider a star formation model dealing with turbulent energy of ISM and an 

Q , averaged approach to take into account the structure of ISM in model of dissipation of turbulent energy. 
^ ■ The proposed models are able to explain the observable star formation pause in our Galaxy corresponding 

' to star ages from 8 — 9 to 10 — 12 Gyr. 

^ ; Star formation model 

^T) ' The main factor which affects star formation is a gas density. However, it is clear that star formation can 

^ be suppressed by supernovae explosions by means of turbulent feedback. Thus, the model star formation 

' law must depend on density and on turbulent energy of gas. 

CN ■ Let us adopt the star formation law in Kennicutt form [H|: 



in 
o 

o 



O I v = c* — 



where p is the gas density, rg cx p"^/^ is the free-fall time, is the star formation efficiency (the values are 
, from 0.1 [IE] to 1 [H]). It is commonly believed that the stars form from a dense and cool gas. It is possible 
^ ■ to adopt this assumption in a simple manner. Just let the star formation efficiency be the fraction 
' of perturbations of density which are unstable by Jeans criterion, assuming the masses of perturbations 
5^ I distributed in a power law: 

c* cx mj^^ , (2) 

where mj cx p"^/^ T^/^ is the Jeans mass [7] depending on density and "temperature" T, and /? is the exponent 
of power distribution of mass of perturbations. It is necessary to explain the meaning of temperature in this 
expression. The temperature value averaged over all components of ISM ~ 10*^K while the temperature of the 
gas in star forming regions is tens and units of Kelvin. However, we can assume the local virial equilibrium 
between thermal and turbulent energy of gas, and the power law dependence of turbulent velocity dispersion 
on scale, so the value of the averaged turbulent energy will reflect the value of temperature and turbulent 
energy on small scales. 

After all substitutions the star formation law will be 

SFR = V' = ffp'^+'rt(i-^) , (3) 
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or, assuming the Salpeter-like exponent P = 2.35, 

2.175 2 

= 5 j^^ME ~ 5 ^2 ' (4) 

where g is the normalizing constant. It may be determined using the star formation law of Tutukov |19] . 
According to this law the star formation process is regulated by ionization and depends on gas density as 

',p = fp^ ^ f = 2xW cm^ g"^ s"^ . (5) 

This law has clear physical base and it works good in the single-zone model of evolution of galaxies [20]. 
The average temperature of modern ISM was chosen to be 4 x 10^ K (though it's value does not affects 
significantly on the chemical abundance history of the Galaxy in single-zone evolution model, see below), 
so the normalizing constant becomes g = 3.2 x 10^^ cm^ g~^ s~^. The temperature T used in the last 
formulae supposed to relate to maximal turbulent velocity dispersion ctq as ctq = k-QT/p,, where fee is the 
Boltzmann constant and /i is the mean molecular weight. 

Some star formation models with turbulent energy account had been offered earlier [HE], but the model 
proposed in this paper has the advantage. This model based on Jeans instability creterion so it can be 
extended for the case when the galaxy rotation, magnetic field or chemical composition is significant. 



Dissipation model 

In all numerical models of galaxy evolution, the ISM in a computational cell presented as a solid 'brick' 
of gas without any structure. However, the ISM has complex turbulent behaviors which are needed to be 
allowed. In a simple approach the ISM may be considered as a set of colliding clouds with different sizes and 
masses. The cloud mass probability distribution function (PDF) 



P{dM} 
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(6) 



and density dependence on scale 



Pi = Po 



1.1 



(7) 



are observable [18] and velocity-scale relation had examined by both observations [12] and computations 

[alE]: 
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Here the scale / is in range {Imin, ^max) and values po and do are averages in an area of size Zmax- 

Imagine that dissipation of turbulent energy occurs trough the collisions of clouds, subsequent contraction 
and heating by shock waves and radiation. It is obvious that the efficiency of dissipation will determine by 
relation of collision time to cooling time [15]. The value of energy radiated away per unit time in this process 
is 



dissipation 
rate 



collision 
frequency 



energy of clouds 
of size I 



fraction of 
radiated energy 



(9) 



The collision rate can be estimated as 
clouds of sizes from / to / -|- dHs 



V3/(2^Gpo 



' P{d/} = pof P{d/} 



Po- 



[T7] . The volume density of turbulent energy of 

(10) 



where vf = ctq — af is the velocity dispersion of clouds of size /. P{d/} is the cloud size PDF which can be 
obtained using the simple assumption about cloud mass and density M; = pil^: 



d/ 
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A = (a - 1)(3 -r) + 1 ^ 1.95 



fll) 
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It is reasonable to describe the fraction of radiated energy for clouds of size / by Poisson distribution 1 — e 
where qi is the relation of the collision time to cooling time with cooling function A: 

qi = TcoiLi 2Tr~ ■ y^^) 

Here the cloud density psh,/ behind a shock wave and the time of collision Tcqu^i can be determined for the 
simple case of head-on collision using the one dimensional approximation, i.e. Riemann problem [11]. In this 
case psh,i = £,iPu where ^; depends on (Tq, Po and and Tcoh,; = //A, where Di is a shock velocity which 

£(1^7—^ where is the dissipation 



is proportional to cloud velocity: Di = ijivi. Finally we can write Q 
efficiency which determines by a complex expression 
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Strong maximum of P{d/} in / 
simplified expression 



Imin, the condition 1^ 



(13) 



1 — exp 



P^ 



Po 



VI VWV^inal)^ 

and the behaviours of ^; and Di lead to 

(14) 



where the single unknown parameter is the maximal scale of turbulence /max- It is set to 5 pc (this is roughly 
the scale of open clusters formation, the turbulence scale cannot be greater than this). 



Star formation in our Galaxy 

The star formation history in galaxies can have a strongly non-monotonic nature showing a couple of 
bursts during a galaxy life time. Those bursts can be stimulated by an accretion of gas from intergalactic 
medium, collisions or close passing of galaxies. The reason to star formation to stop can be the supernovae 
explosions P]. In star formation history of our Galaxy an epoch exists when the star formation process 
was suppressed. It is seen with distribution of [Fe/0] fB], [Eu/Ba] - [Fe/H] [T3] and [Mg/Fe] - [Fe/H] [5J. 
This epoch may be interpreted as a pause between the end of formation of thick disk and the beginning of 
formation of thin disk [SJ [T3|. The pause corresponds to star ages between 8 — 9 and 10 — 12 Gyr. 

The model of star formation and of dissipation proposed above, was taken to explain this pause in a 
hierarchical scenario. These models was built into the single-zone model of evolution of galaxies (Firmani & 
Tutukov [4j) which was initially implemented by Wiebe [201 • The detailed description of modified model and 
target setting has given in [TO]. In modified model the star formation suppressed by collision of Galaxy and a 
satellite of mass 2 x 10^'^ Mq and radius 6.32 kpc at redshift z = 1.5. The results of modelling are presented 
on Fig. [H After collision the temperature increases almost by two orders causing decrease of dissipation 
efficiency (see with the result that star formation process ceases due to high temperature (see @) for 
a 1.5 Gyr. The pause in star formation process is seen as a plateau on the abundance graphs [O/H] and 
[Fe/H] (solid line). 

It is interesting to compare the abundance history with distribution of star abundance in Solar neighbor- 
hood. On Fig. (I2J the [O/H] — [Fe/0] distribution is shown for near metal-poor stars f6], and the abundance 
history in hierarchical scenario. The arrow marks the pause of star formation which is close to the bound of 
populations of thick disk and thin disk. The consequence of this pause is clear: we can see the discontiniuty 
of the stars age with continuous abundance history (see the plateau on Fig. [T|). The standard model of 
Firmani & Tutukov [4| does not give such a plateau, giving the right position and shape of the evolution 
track, though. 



Conclusions 

The need of taking into account the turbulent energy in star formation model is obvious. The large differ- 
ence of properties of ISM on different scales is obvious also. In the present paper the model of star formation 
and dissipation of turbulent energy was offered and implemented to single-zone model of galaxy evolution. 
It is shown that the star formation history and history of chemical abundance in Solar neighborhood can be 
modeled in the frame of this model using the scenario of collision of the Galaxy and the small satellite. 
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Fi gure 1: The results of modelling of collision of Galaxy and small satellite. SFR is the star formation rate in solar masses per year; T 
is the temperature; is the efficiency of dissipation. Dotted line represents the standard model of Firmani & Tutukov [1], dashed line 
traces the evolution of a satellite before collision and solid line shows the evolution of the collision area in Galaxy, after collision. 
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Figure 2: Evolution of chemical abundance of Galaxy in hierarchical scenario. Dot-dashed line is the standard model of Firmani & 
Tutukov d], solid line is the satellite and dashed line is the collision area in Galaxy. The arrow marks the position of pause in star 
formation process which is seen as a plateau on the abundance graphs on Fig. [T] Symbols are halo stars (open squares), stars of the thick 
disk (filled circles) and thin disk (open triangles). The stars was classified using accurate stellar kinematics [6]. Error bar is at bottom 
right (there are same error bars for all the stars). 
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